####################################################
# 
# --- Command log file ---
# PLAXIS 2D 2025.1
# Tutorial - 16
# Pile driving
# Note: Use the PLAXIS Commands runner to execute this
# 
####################################################
_initializerectangular SoilContour 0 0 30 18
_setproperties "Title" "Pile driving" "Company" "Bentley Systems Inc" "Comments" "" "UnitForce" "kN" "UnitLength" "m" "UnitTime" "day" "UnitTemperature" "K" "UnitEnergy" "kJ" "UnitPower" "kW" "WaterWeight" 10 "ReferenceTemperature" 293.15 "LiquidSpecificHeatCapacity" 4181.3 "LiquidThermalConductivity" 0.0006 "LiquidLatentHeat" 334000 "LiquidThermalExpansion" 0.00021 "LiquidTemperature" 293.15 "IceSpecificHeatCapacity" 2108 "IceThermalConductivity" 0.00222 "IceThermalExpansion" 5E-5 "VapourSpecificHeatCapacity" 1930 "VapourThermalConductivity" 2.5E-5 "VapourSpecificGasConstant" 461.5 "UseTemperatureDepWaterPropsTable" False "ModelType" "Axisymmetry" "ElementType" "15-Noded"
_borehole 0
_soillayer 0
_soillayer 0
_set Soillayer_1.Zones[0].Top 18
_set Soillayer_1.Zones[0].Bottom 7
_set Borehole_1.Head 18
_soilmat "Identification" "Clay" "SoilModel" "Mohr-Coulomb" "DrainageType" "Undrained B" "Colour" 15262369 "gammaUnsat" 16 "gammaSat" 18 "ERef" 5000 "nu" 0.3 "EInc" 1000 "verticalRef" 18 "sURef" 5 "sUInc" 3 "InterfaceStrengthDetermination" "Manual" "Rinter" 0.5  "K0Determination" "Automatic"  
_soilmat "Identification" "Sand" "SoilModel" "HS small" "DrainageType" "Undrained A" "Colour" 10676870 "gammaUnsat" 17 "gammaSat" 20 "E50Ref" 50000 "EOedRef" 50000 "EURRef" 150000 "G0Ref" 120000 "gamma07" 0.0001 "phi" 31  "UseDefaults" True "K0Determination" "Automatic"      
_soilmat "Identification" "Pile" "SoilModel" "Linear Elastic" "DrainageType"  "Non-porous"  "Colour" 10283244 "gammaUnsat" 24 "ERef" 30000000 "nu" 0.1 "InterfaceStrengthDetermination" "Rigid" "K0Determination" "Automatic"    
_setmaterial (Soil_2) Sand
_setmaterial (Soil_1) Clay
_gotostructures
_polygon 0 18 0.2 18 0.2 7 0 7
_delete (Polygon_1)
_line 0.2 18 0.2 7
_undo
_line 0.2 6.6 0.2 18
_neginterface (Line_1)
_lineload 0 18 0.2 18
_set Line_2.LineLoad.LineLoad.qy_start -1
_loadmultiplier
_set Line_2.LineLoad.LineLoad.Multipliery LoadMultiplier_1
_set LoadMultiplier_1.Amplitude 5000
_set LoadMultiplier_1.Frequency 50
_mergeequivalents Geometry
_gotomesh
_mesh 0.06
_gotostages
_phase InitialPhase
_set Model.CurrentPhase Phase_1
_set Phase_1.Identification "Pile activation"
_activate (NegativeInterface_1_1) Phase_1
_set Soil_1_1.Material Phase_1 Pile
_phase Phase_1
_set Model.CurrentPhase Phase_2
_set Phase_2.DeformCalcType "Dynamic"
_set Phase_2.Identification "Pile driving"
_set Phase_2.Deform.TimeIntervalSeconds 0.01
_set Phase_2.Deform.ResetDisplacementsToZero True
_set Phase_2.MaxStepsStored 100
_activate (DynLineLoad_1_1) Phase_2
_set Dynamics.BoundaryXMin Phase_2 "None"
_set Dynamics.BoundaryYMin Phase_2 "Viscous"
_phase Phase_2
_set Model.CurrentPhase Phase_3
_set Phase_3.Identification "Fading"
_set Phase_3.DeformCalcType "Dynamic"
_set Phase_3.Deform.TimeIntervalSeconds 0.19
_set Phase_3.MaxStepsStored 100
_deactivate (DynLineLoad_1_1) Phase_3
_selectmeshpoints
/output addcurvepoint "Node" Soil_1_1 (0 18)
/output update